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Abstract. We consider a derivative nonlinear Schrodinger equa- 
tion with a general nonlinearity. This equation has a two param- 
eter family of solitary wave solutions. We prove orbital stabil- 
ity/instability results that depend on the strength of the nonlin- 
earity and, in some instances, their velocity. We illustrate these 
results with numerical simulations. 

1. Introduction 
The derivative nonlinear Schrodinger (DNLS) equation 
(1.1) id t u + d 2 x u + i(\u\ 2 u) x = 0, 

is a nonlinear dispersive wave equation that appears in the long wave- 



length approximation of Alfven waves propagating in a plasma 25 , 26 
30] . Applying the gauge transformation 



^1.2) < ' = */(./•) cxp / ^ - j \ii{ii)\ 2 (hj 



equation (1.1) has the form 

(1.3) id t tp + d 2 x tp + i\tfj\ 2 ^j x = 0, xeR. 

This equation has a Hamiltonian structure and can be written as 

dtp 



at ^ 



where the Hamiltonian is 



1 f°° 1 
E = -j \iP x \ 2 dx + -% j | ( frr,,/.r. 
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Equation (1.3) and some of its generalizations, also appear in the mod- 
eling of ultrashort optical pulses (T 27 . Furthermore, the DNLS equa- 
tion has the remarkable property of being integrable by inverse scat- 
tering [18] . It admits a two-parameter family of solitary wave solutions 
of the form: 
(1.4) 

c, ,3 
4 



u u , c (x, t) = tfiu jC (x - ct) exp % \ out + -(x - ct) 



x—ct 



<pl,c 



where u > c 2 /4 and 
(1.5) <PuM = 

is the positive solution to 




(4w - c 2 ) 



h(a\/Auj — c 2 y) 



C N 



1.6) 



+ (u 



16 



0. 



Guo and Wu 11 showed that these solitary waves are orbitally stable 
if c < and c < Au. Colin and Ohta |2| subsequently extended the 
result, proving orbital stability for all c, cr < 4a;. 



Definition 1.1. Let u W;C be the solitary wave solution of (1.1). The 
solitary wave u u>c is orbitally stable if, for all e > 0, there exists 5 > 
such that if ||wo — UujAIh 1 < 8, then the solution u(t) of (1.1) with 
initial data u(0) = Uq, exists globally in time and satisfies 



sup inf 

t>0 (0,2/)6Txl 



\u(t) -e i0 u u , c (t,- -y)\\m < e. 



Otherwise, it WiC is said to be orbitally unstable. 

In an effort to understand the structural properties of DNLS, we 
study an extension of (1.3) with general power nonlinearity (cr > 0). 

(i.7) i^v* + ^ + = o. 



Equation (1.7) also admits a two-parameter family of solitary wave 

solutions, 

(1.8) 



c 1 

Vvfo *) = w(z-ci) exp i <J ut + - (x - ct) - 2a + 2 



x—ct 



9Vc 



(r/)rfr/ 



where a; > c 2 /4 and 
(1.9) ^, c (y) 2CT 



( ( x+l)(4o;-c 2 ) 



2 v ^(cosh(crV4o; - c 2 y) - ^) 
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is the positive solution of 
(1.10) 

d2 , / ^ \ , I l2rr "I" 1 l l/lrr n 

y ¥W + (w- j )<p u ,c+ glwl <Puj,c- r 2 a + 2)^ (fw ^ Vu,c = 

It is convenient to define 

(1.11) = ^c{yy 9uiAy \ 

with the traveling phase 

(1.12) e u M = c - y - -1— [ v ^Mdv- 

2 2a + 2J_ 00 

Clearly, 

(1.13) < c (a:,*) = e <wt ^, c (x-cf) > 
and the complex function <p u , c (y) satisfies 

(1.14) - dy^^ + u<j) UiC + idd x 4> u ,o ~ «|<?W| 2<T <9<,0^ C = 0,1/6 1. 

Provided there is no ambiguity, we write 0, for WjC , <£> WjC respec- 
tively. Furthermore, we only consider admissible values of (u, c) satis- 

2 

fying the conditions oo > x> c6l. 

1.1. Main Results. We investigate the stability of solitary wave solu- 
tions ^ WiC to the gDNLS equation. This is determined by both the value 
of a and the choice of the soliton parameters, c and oo. These results 
are conditional in the sense that for a ^ 1, we lack a suitable local well- 
posedness theory. Throughout our study, we assume that given a > 0, 
and ip G H l (R), there exists a weak solution ^ eC ([0,T); H l (R)) of 
(1.7), for T > 0, which satisfies 

(1-15) f t (^(;t),f) = (E'^(.,t)),-Jf) 

for appropriate test functions /. E is the energy functional, and J is 
the symplectic operator. These are defined in Section [2] 
Subject to this assumption, we have the following results: 

Theorem 1.2. For any admissible (oo,c) and a > 2, the solitary wave 
solution i[> u> c{x,t) of (1.7) is orbitally unstable. 

For a between 1 and 2, slow solitons, those with sufficiently low c, 
will be stable while fast right-moving solitons will be unstable: 

Theorem 1.3 (Numerical). For a G (1,2), there exists z = Zo(a) G 
(—1,1) such that: 
(i) the solitary wave solution %l) Wfi {x,t) of (1.7) is orbitally stable for 
admissible (lj, c) satisfying c < 2z y/u. 
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(ii) the solitary wave solution if)u} tC (x,t) of (1.7) is orbitally unstable 
for admissible (lj, c) satisfying c > 2z Q y/uj. 



Our last result concerns a < 1, where all solitons are stable: 

Theorem 1.4 (Numerical). For admissible (u, c) and < a < 1, the 

solitary wave solution ^^(a^t) of (1.7) is orbitally stable. 



The endpoint, a = 1, corresponds to the cubic case, which has al- 



ready been studied in [2 


11 . 


Theorems 


1.3 


and 


1.4 


are fully rigorous up to the 



of the sign of a function of one variable, which is parametrized by 
a. The number z in Theorem 1.3 corresponds to a zero crossing. 
This function, defined below by (4.3), includes improper integrals of 
transcendental functions. 



Theorem |1.3| is particularly noteworthy for distinguishing gNLS from 
the focusing nonlinear Schrodinger equation (NLS), 

(1.16) 



iipt + A</>+ |^| 2ff ^ = 0, V 



Equation (1.7) is invariant under the scaling transformation ijj\(x,t) = 
'ip (Xx, X 2 t) . This implies that its critical Sobolev exponent is s c = 
it is L 2 -critical for a = 1, and it is L 2 -supercritical, 



\ — 77-. Hence, 

energy subcritical for cr > 1. While NLS only admits stable solitons 
in the L 2 -subcritical regime, gDNLS admits stable solitons not only in 
the critical regime, but also in the supercritical one. 



1.2. Outline. Our results are proven using the abstract functional 
analysis framework of Grillakis, Shatah and Strauss 



sec 



34 35 



for related results and 31 for a survey. The test for stability involves 
two parts: (i) Counting the number of negative eigenvalues of the lin- 
earized evolution operator near the solitary solution of (1.7), de- 
noted n(iJ^); (ii) Counting the number of positive eigenvalues of the 
Hessian of the scalar function d(u, c) built out of the action functional 
evaluated at the soliton, denoted p(d"). We give explicit characteriza- 
tions of and d(cu, c) in Section |2] We then apply: 



Theorem 1.5 (Grillakis et al. [7|[8]). 

(1.17) p(d") < n{Ht) 

Furthermore, under the condition that d is non- degenerate at (u,c): 
(i) If p(d") = nlHfj,), the solitary wave is orbitally stable; 
(ii) If n^H^) — p(d") is odd, the solitary wave is orbitally unstable. 

In Section |3j we show that n^H^) = 1, and in Section [ij we compute 
the number of positive eigenvalues of the Hessian, d". We complete the 
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proofs of Theorems L2 to 1A in Section |5j In Section [6j we discuss the 
results and illustrate them with numerical simulations of both stable 
and unstable solitons. Algebraic manipulations useful for the evalua- 
tion of det[d"(u>, c)] and tr[d"(co, c)] are presented in the Appendix. 



1.3. Remarks on Well Posedness Assumption. The DNLS equa- 
tion (cubic nonlinearity) has been studied in H l and higher regular- 
ity spaces, with results for both local and global well posedness re- 
sults 



13 16,29 



tion related to 



32 



33 



Much of the analysis relies on a transforma- 
\.2\ that turns the equation into two coupled semi- 

In addition, Hayashi 



linear Schrodinger equations with no derivative 
T3lfl4l[l6l identified a smallness condition on the data 



'1.18) 



\Uq £2 < V27T 



for which global solutions exist in H s , seN. The constant y27r is the 
L 2 -norm of the ground state of the quintic NLS soliton. More recently, 



the global in time result for data satisfying (1.18) were extended to H s 
spaces with s > 1/2 in (3). DNLS with low regularity has also been 
studied on the torus, (9l[l0l. 



There has also been progress beyond the cubic equation in the afore- 
mentioned results. Some studies, such as (6,29 32 , include additive 
terms to the cubic nonlinearity with derivative. More generally, Kenig, 



Ponce and Vega 19,20 used viscosity methods to show that, for gen- 



eral quasilinear Schrodinger with polynomial nonlinearities, local well- 
posedness holds in Sobolev spaces of high enough index (See Linares- 
Ponce [22] for a review). In |12|, Hao proved that ( 1.3[ ) is locally 
well-posed in H l l 2 intersected with an appropriate Strichartz space 
for a > 5/2. Working in the Schwartz space, Lee [2l] used the frame- 
work of inverse scattering to show that DNLS is globally well-posed for 
a dense subset of initial conditions, excluding certain non generic ones. 



2. Problem Setup 

In this section, we define the linearized operator, the invariant quan- 
tities and the action functional of the soliton. Throughout, we adopt 
the notation of [8]. 

We study the problem in the space X = if 1 (IR), with real inner prod- 
uct 

(2.1) (u, v) = ^R. / (u x v x + uv)dx. 
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The dual of X is X* = H~ 1 (R). Let / : X -> X* be the natural 
isomorphism defined by 

(2.2) (Iu,v) = (u,v), 
where (•, •) is the pairing between X and X*, 

(2.3) (f,u) = 3? / fudx. 

Jn 

gDNLS can be formulated as the Hamiltonian system 

(2.4) f = JE> W , 

where the map J : X* — > X is J = — z, and .E is the Hamiltonian: 

(2.5) E=- \ip x \ 2 dx+ 3 / iVf^^Mx. (Energy) 
Two other conserved quantities are: 

(2.6) Q = - / |V| 2 rfx, (Mass) 



oo 
oc 



1 /"°° - 

(2.7) P = — S / il)ip x dx. (Momentum) 

2 J -oo 

Let Ti and T 2 be the one-parameter groups of unitary operator on X 
defined by 

(2.8a) Ti(t)$(x) = e-' l %(x), $(x) E X, i G 1 

(2.8b) T 2 (t)$(x) = $(x + t). 



Then 
and 



ipw,c = Ti(-ujt)T 2 (-ct)(f) U)C (x). 



9 t Ti(wt)| t=0 = -iu, d t T 2 (ct)\ t=0 = cd : 
We define the linear operators 

(2.9a) B 1 = J- 1 d t T 1 (ojt)\ 1 ^ = u, 

(2.9b) B 2 = J- 1 ^^)!^ = i^. 



The mass and momentum invariants (2.6) and (2.7) are thus related to 
the symmetry groups via: 

(2.10a) Qi ee ^(B 1( f>, 4>) = \®J uWdx = | J |0| 2 cfe = coQ, 
(2.10b) Q 2 ee ^(B 2 4>,4>) = J icdxWdx = -ic3 J d x ^dx = cP. 
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Computing the first variations of (2.5), (2.6) and (2.7), we have 

(2.11a) £'(0) = -<9^-i|0| 2CT d^ 

(2.11b) Q'(0)=0, P'{4>) = idA- 

The second variations are: 



(2.12a) 
(2.12b) 



E"{<t>)v 



-d x - icr\ 



I 2ct- 2" 



%<P-t\<P\ 2tT d x )v 



1 2cr-2. 



Q"((j))v = v, P"((j))v = id x v. 



2.1. Linearized Hamiltonian. The linearized Hamiltonian about the 
soliton d> is 



(2.13) 



H <j> u = [E"(cf>)+Q , ((4>)+Q 2 , (4>)}u 

= [E"{<P)+ujQ"{<P) + cP"{<P)}u 



for u e H 2 (R). For later use, we give two equivalent expressions of H^. 
First, we decompose it into complex conjugates: 

Lemma 2.1. For any function u in the domain of H^, 

H^u = Liu + L 2 u, 

where 



(2.14a) 
(2.14b) 



L\ = —d 2 . + to + icd x — ia\(j>\ 



2cr-2, 



I 2(7 



0,- 



L 2 = —ia\ 



|2(7-2 



Second, we give the expression of and the quadratic form it in- 
duces, after extraction of the soliton's phase: 

Lemma 2.2. Let u 6 H 1 be decomposed as 



(2.15) 



u = e l6 (ui + iu 2 ) 



where 9 is given by (1.12), and U\ and u 2 ore the real and imaginary 
parts of ue~ %e . Then 



(2.16) 
and 

(2.17) 



H^u = e %e [(Luui + L 2 iu 2 ) + i {L u ui + L 22 u 2 )] 

{H^u.u) = (Luu^ux) + (L 2l u 2 ,ui) 

+ (L12UUU2) + (L 22 u 2 ,u 2 ), 
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c 2 c(2a + 1) 4a 2 + 6a + 1 
(2.18a) Ln^-^ + c-^ + ^-^y 2 -- 



(2.18b) L 21 
(2.18c) L 12 



4(7 



a 



,„2(7-l,„ i ° ,„2a£> 



(2o-+l)<7 2a _ x a 



(7+1 



a + 1 



(2.18d) L 22 ^-d yy + u- C - + l^-^±^^. 



Proof. Using ( |1.12[ ) 
(2.19) 



0,, 



2 2a + 2 



2(7 



a 



2(7-1 



a + 1 



We can then rewrite the operators L\ and L 2 from Lemma 2+ in terms 
of <p as 



Li = -dy +u + zc9y - 2y? 2<T <9 ; 



CTC 2a CT 4<7 



2 Y 2a + 2 



<p™ - zay2 2<T 1 (p y , 



CO 



-if 



:2a 



a 



4(7 



,2(7-1, 



2 r 2a + 2 
Letting x = ue~ l °, we have 
J0 



„2i0 



H^u =e v 



8 +u ° 2 I C ^ +1 ) ^ ^ ..2(7-1 



(2.20) 



/a 



a+l^ l(a+l)^ 



2a 2 + 4a + l 4(T 

V 5 



X 



+ e Ji 



c:a 



-<p 



2a 



a 



2(7-1, 



2a + 2 



ip — tcrip (p y 



Since \ — u i + ^ M 2, we get (2.16) and (2.17) by grouping terms appro 
priately. 



□ 



2.2. Scalar Soliton Function. Using (1.14), we observe that when 
evaluated at the soliton 0, 

(2.21) E' + cuQ' + cP' = 0. 
For any u > c 2 /4, we define the scalar function 

(2.22) d{u, c) = E(<f> u , e ) + Qi{<f> u , c ) + <3 2 (<M, 

which is the action functional evaluated at the soliton. It has the 
following properties: 
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Lemma 2.3. 

(2.23) 
(2.24) 
(2.25) 

The Hessian is 
(2.26) 



d(uj, c) 
dud(u, c) 
d c d(uj, c) 

d"{u,c) 



£(</>)+ wQ(0) + cP(0), 
Q{<t>) > 0, 

P{4>). 



d^p^) d c P(<f>) 



Proo f. Us ing ( |2.10a[ ), ( |2.10b[ ) and ( |2.22[ ), we have ( |2.23[ ). Di fferen tiat- 
ing (2.23) with respect to u and c respectively and using (2.21), we 
obtain (2.24) and (2.25). The expression for the Hessian follows. 

□ 



3. Spectral Decomposition of the Linearized Operator 

This section provides a full description of the spectrum of the lin- 
earized operator H^. In particular, we prove: 

Theorem 3.1. For all values of a > and admissible (lu, c), the space 
X = H 1 can be decomposed as the direct sum 

(3.1) X = N + Z + P, 

where the three subspaces intersect trivially and: 

(i) N is a one dimensional subspace such that for u £ N , u ^ 0, 

(3.2) (H^u,u)<0. 

(ii) Z is the two dimensional kernel of H^. 
(Hi) P is a subspace such that for p G P, 

(3.3) (H^p,p)>5\\p\\ 2 x 
where the constant 5 > is independent of p. 

Corollary 3.2. For all values of a > and admissible (oo,c), 

n(H^) = 1. 

An important ingredient of the proof involves rewriting the quadratic 
form (2.17) induced by in a more favorable form. This rearrange- 
ment, inspired by [TT] , expresses it as a sum of a quadratic form involv- 
ing an operator with exactly one negative eigenvalue and a nonnegative 
term. 
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Lemma 3.3. Let 



u = e l9 (ui + iu 2 ), 



where 9,Ui,u 2 are defined the same as Corollary \2.2\ then 
(3.4) 



{H^u.u) = (Lxxu\,ui) 
where 



P ((p l u 2 ) + 



a 



(7+1) 

8a 2 + 6a + 1 



dy, 



la 



(3.5) 

Proof. Recall the terms in the quadratic form ( 2.18[ ). We first examine 
L\\. The relationship between Lu and L n is 

■2 



(3.6) 



L 



ii 



L n + 



a 



+ 



rP 



la 



Next, consider L 22 . From (1.10), 
(3.7) L 22 ip = 0. 

Letting u 2 = (p~ l u 2) we can then write 
(3.8) 

/ c 2 c 

(L 22 u 2 , u 2 ) = (-d yy u 2 , u 2 ) + ((uj-j + -tp™ - + 



2a 2ff + 1 4aN 

<P )u 2 ,u 2 



[-<PyyU 2 - 2(PyU 2 y ~ fU 2 yy , <£U 2 ) 



// C C 

+ ((^-^ + 2^ 



2a 



2(7 + 1 



■V9 4<T )V3M2, <PU 2 ) 



4(a + l) 2 

= (U 2 L 22 lf, (fU 2 ) + (-2(f y U 2y - <fU 2 yy, (fU 2 ) 
= (-(ip 2 U 2y )y,U 2 ) = (<pU 2 y,<pU 2 y}, 

where u 2y and u 2yy denote d y u 2 and d yy u 2 , respectively. Lastly, we 
simplify the off diagonal entries, L 2 i an d L 12 . Integrating by parts, we 
have 



(L 12 ux,u 2 ) 



(2(7 + l)(7 



,2a-l 



(7 + 1 
2(7 + 1 



— f7 —-(p 2a d y )u 1 ,u 2 
a + 1 y J 



<(^ 2<T ) y ,M 1 M 2 ) - — ^-j- (<f 2a U ly ,U 2 ) 



2((7 + l) 

2(7 + 1 , 2 > 2(7 + 1 , 2 \ (7 , 2 

{If U 2 y, Ul ) + ^^^^ {<P U ly , U 2 ) - {(f U ly . 



2((7+l) ^ ^ L/ 2((7+l) 
= 27aTl) <^ W2 - Wl > + 2(a+T) <^«"«*> ' 



(7 
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Similarly, 

(L 2 iM 2 , «l) 



" 2( g + i) < ( »' 2 ° , »' Kl " 2 ) + ^TT (^""^"i) 

1 1 (T 

2( ff + i) <^ 2,J ^' Ul > + 2 (g + i) ^ 2CTui ^ U2 > + ^TT <^ 2,J ^' Ul > 



1 



The off diagonal terms then sum to 
IT \ It v 2(7 + 1 



1 



— Up 2<T U2y, u x ) H — Up 2<T u Xy , M2) 

1 cr + 1 



(T+ 1 ' (re- 

introducing £t2 = v 9_lu 2 into the above expression, and integrating by 
parts, 

2a 



(3.9) 



{L l2 ux,u 2 } + (L 2 iu 2 ,ui) 



a + 1 



ip 2a+1 U 2 y,U 1 ) 



Combining (3.6), (3.8) and (3.9) 



{H^u.u) =(LnU 1 ,«i) + (<pu 2y ,(pu 2y ) 



a 



:(LnM!,Mi) 



(7 



-2(7 



a + 1 



V ui ) + 



<PU 2 y + 



a 



2a 

a + 1 

2 



<p 2C Ul,tpU 2 y 



a+1 



dy. 



□ 



3.1. The Negative Subspace. Next, we characterize the negative 
subspace, N. For that, we need the following lemma on L n . 

Lemma 3.4. The spectrum of L\\ can be characterized as follows: 

• L\\ has exactly one negative eigenvalue, denoted — A 2 1; with 
multiplicity one, and eigenfunction xn, 

• G a(Lu), and the kernel is spanned by ip y , 

• There exists fin > such that 

a{L 11 ) \ {-A^O} C [/in, 00). 

Proof. First, we observe that since p is exponentially localized, Ln is 
a relatively compact perturbation of — d y + u — \. By Weyl's theorem, 
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the essential spectrum is then 



u - ^,00 



Consequently, all eigenvalues below the lower bound of the essential 
spectrum correspond to isolated eigenvalues of finite multiplicity. By 
differentiating (1.10) with respect to y, we see that 

(3.10) Ln(p v = 0. 

Hence, L n has a kernel. Viewed as a linear second order ordinary 
differential equation, Lnf = has two linearly independent solutions. 
As y — > —00, one solution decays exponentially while the other grows 
exponentially. Thus, up to a multiplicative constant, there can be at 
most one spatially localized solution to Lnf = 0. Therefore, the kernel 
is spanned by <p y . 

From Sturm-Liouville theory, this implies that zero is the second 
eigenvalue of Ln, and Ln has exactly one strictly negative eigenvalue, 
— A n , with a L 2 normalized eigenfunction xii- 



(3.11) £11X11 = -^liX 



11 • 



If we now let 



(3.12) fin = inf 



£11/,/ 



f^0J± Vy J±xn (/, /) 

we see that fin > 0, since if it were not, it would correspond to another 
discrete eigenvalue less than or equal to zero. It is either a discrete 
eigenvalue in the gap (0, u — x) or the base of the essential spectrum. 
Regardless, cr(Ln) \ {— A n , 0} is bounded away from zero. 



□ 



Using xn, we construct the negative subspace N. 

Proposition 3.5. Let 

(3.13) iV = span{x-} 

where 

(3.14a) x- = (Xn + % 2 )e ie , 

/y 
^ 2a - 1 (s) X n(s)ds + k 12 
_ 

and ki2 G 1 is chosen such that 

(3.15) (Xi2,^) = 0. 
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(H^,u,u) < 0. 
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Proof. The function X12 is in L 2 . Indeed, the integral in (3.14b) is well 
defined since, as \y\ — > oo, 

\<p(y)\ <exp|-v/o;-c 2 /4|2/|} , 

\xu{y)\ ^ ex P |-\/w-c 2 / 4 + ^i Iz/lJ ■ 

Thus the integrand is bounded. 
From (|Ofc and (J3TTT1), 



(^X-,X-) = (£11X11, Xn) = -An < 0. 



□ 



3.2. The Kernel. In this subsection, we give an explicit characteriza- 
tion of the kernel of H^. 

Proposition 3.6. Let 

(3.16) Z = span{xi,X2} 

where 



(3.17a) 



Xi=[<Pv + ^2 ~ o — a )<P e 



2a + 2 



(3.17b) X 2 = i^e id 

with k 2 is a real constant such that 



(3.18) 

Then Z = ker Ha. 



<(F)<p,<p)=0. 



2a + 2' 



Proof. We first prove that xi and X2 are linearly independent elements 
of the kernel, and then show that the kernel is at most two dimensional. 
Applying (in the form (2.16)) to X2 and using that L 2 itp = and 
(3.7), we get H^x? = 0. For xii we compute 

L n <f y + L 21 (k 2 - 0j 1 n <p 2lT )<p 



2a + 2' 



2ct+1 



a~ 



.la 



(a + iy 



1 2a 2 
2a + 2a + l 



if^iPy = 
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and 



L 12 ip y + L 22 {k 2 - - f 2c7 )f 
2(7 + 2 

2a 2 + a , a 



,Jla-\,Jl w ,„2cr 

V fy ~T<P Vyy 



(7 + 1 r ry (7 + 1 

' [-2(7(2(7 + l)^ -1 ^ - 2<VV ra + ^ 2ct L 22 ^] = 0. 



2a + 2 

Thus, Z C kerif^, and the kernel is at least two dimensional. 

We now show that it is exactly two dimensional. If we consider the 
problem 

H+f = 0, 

as a second order system of two real valued functions, we know there are 
four linearly independent solutions. As y — > — oo, two of these solutions 
decay exponentially, while two grow exponentially. Thus, there are at 
most two linearly independent solutions which are spatially localized. 
Hence, Z = ker H^. 

□ 

3.3. The Positive Subspace and Proof of the Spectral Decom- 
position. We define the subspace P and prove Theorem 3.1 For that, 
we need the following lemmas about L\\ and L 22 . 

Lemma 3.7. For any real function f G if 1 (IR) satisfying the orthogo- 
nality conditions 

(3.19) (/,¥>„> = (/,Xn) = 0, 
there exists a positive number <5 n > 0, such that 

(3.20) (W,/)>M/II^- 



Proof. From Lemma 3.4, (3.12) holds on the subspace orthogonal to 

ip y and xu, so 

(w,/)>mimiL- 

To get the H 1 lower bound, let 

x r / \ c 2 c(2(x + l) 2CT 8(t 2 + 6(t + 1 4<t 

^ = w " 4 + - 4((7 + l) 2 * > 

so that Ln = —d yy + Vi, with || Vi || j^oo < oo. Thus, 

(L u f, /} = (-d yy f, /) + (V x f, f) 

><-^y/./>-||Vi||L-||/||Li 

= R/il^- — INMW,/}- 
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It follows that 



(L n f,f)> 



WW) 



1 + tin || Vi\\ L <x> 
Taking 5 U sufficiently small, we have 



(L n f,f)>5 n 



H 1 ' 



□ 



Lemma 3.8. For any real function f 6 H 1 ^) satisfying 

(3.21) l/» = 0, 

t/iere exists a positive number 622 > 0, such that 

(3.22) (lW,/)>M/||iri. 

Proof. As was the case for L n , L 2 2 is a relatively compact perturbation 
of — dy + co — c 2 /4, so it also has 



&CSS 



22 I 



OJ — ^, OO I . 



Thus, all points in the spectrum below u — c 2 /4 correspond to discrete 



eigenvalues. From (3.7) and if is strictly positive, Sturm-Liouville the- 



ory tells us that zero is the lowest eigenvalue. Let 

(3.23) /i 22 = mf . . 

(/,/) 

We know that /i 2 2 > 0, otherwise this would contradict with being 
the lowest eigenvalue. Therefore 

(w,/>>/i 22 imi L2 . 

Using the same argument as in Lemma 3.7, we obtain (3.22). □ 



We now prove Theorem 3.1 



Proof. Recall N, Z as defined define by (3.13) and (3.16). We define P 

as 

(3.24) 

P = {peX I 0R(e-^) )X ii> = (n(e- id p), Vy ) = (9f( e -»p), ^) = 0}. 
We express u G X as 

u = axx- + (bxXi + &2X2) + p, 

where 



«i = (ui,Xu), b i 



2 ' 



\<Pv\\h 



(v<2,<P) 
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with Mi and u 2 are real and imaginary part of e~ ld u. Clearly, a\x~ G N 
and 61X1 + &2X2 G Z. It suffices to show p G P. We write p = 
{pi + ip2)e ld with pi and P2 real. Since (p y is odd and xn is even, 
(fyiXii) — 0, and we readily check that (pi,Xn) = (pi^y) = 0. 
Furthermore, by (3.15) and (3.18), we also have (p2,<f) = 0. Thus, 
p6P, and u is indeed decomposed into elements of N, Z and P. 
Finally, we show that is positive on P. Let p2 = ip ~ 1 P2- By (3.4), 

a 



(3.25) {H^p.p) = (L n pi,Pi) 



((pdyp 2 



a + l 



Lemma 3.7 gives the desired lower bound on the first term. For the 
second term, we break it into two cases, depending on how ||v5^p 2 || L2 
and ||pi|| L 2 compare. Let 

2(7 



(3.26) 



1 2(7 



a+l 



M\L° 



(a) If IIv^^s/^IIl 2 > Co-lbilU 2 ; we estimate the second term in (3.25) as 
follows, 



fdyP2 



a 



a 



a 



L 2 



a + l 



I I|2ct 11 || 
Mli« |bl|| L 2 



\\<fd y p 2 \\ L2 - -c a \\piWv > ^ ydyp 2 \\ L 2 



By (3.8), we then have 



fdyP2 + 



a 



a + l 



^Pi 



L-' 



> J (L 22 P2 1 P2) 



By Lemmas |3.7| and 3.8, we get 



(H4,p,p) > (ZnPuP^ + - (L 2 2P2,P2) > $a \\P\\ 2 H 1 , 

for some small enough 5 a . 
(b) If instead, U^d^lU 2 < Co-|bi||z, 2 , then, 

(H^p) > (ZnPuPi) > -ylbill^ + ^ylbilli 2 



(3.27) 



6l -\\Pi\\m + ^k(L 2 2P2,P2) > 8b \\pf H i ■ 



2 ' 2C G 

Taking the smaller value of 6 a and 5^ as 5, we have 

(3.28) (H (t> p,p)>S\\p\\ 2 H1 . 
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It follows that N, Z and P have trivial intersection amongst one an- 
other. Hence X = N + Z + P. □ 

4. Analysis of the Hessian Matrix 

In this section, we compute the number of the positive eigenvalues 
of the Hessian matrix of d(u, c), p(d"(u), c)). Since the number of neg- 
ative eigenvalues of H$ is in all cases equal to one, p(d"(u,c)) will 
determine whether or not the soliton is stable. 

To make this assessment, we examine the determinant and the trace 
of d"(u,c). From Lemmas 2.3 and A. 3 the determinant can be ex- 
pressed as 

det[d"{u, c)] = d u Qd c P - dcQd^P 

= 2-- - V -2 (l + a)*{4u - c*)°- % cj-«- 2 

x [4(cr — l)ua — 2y/uca + (4a; — c 2 )ai] 

x [4(cr — 1)o;q!o + 2-y/a7cao — (4a; — c 2 )cti] , 



(4.1) 



where 



a n {u), c; cr) 



/r 



> 0, 



cos 



h(o-V4 



a; — crx 



c 

2^' 



cr; Co>, c) 

Meanwhile, the trace is 

ti[d"(u,c)] = d ul Q + d c P 

(4.2) = 2-°- 2 a-\l + tr)£(4w - c 2 )- _1 (l + uj)oj~^ 

x (c(c 2 - 4w)ai + 2 v / w(c 2 - 4(cr - l)u>)a ). 

Theorem 4.1. If a >2, and > c 2 , p(d"(u,c)) = 0. 

Proof. We examine the terms appearing in ( |4.1 ). The first term is 
clearly positive. The second term is also positive, 

4(cr — l)ua — 2^fojca.Q + (4u — c 2 )ai 

c. r 2 

a- 1 



2^ )ao+{1 -h )ai 



=Auj 

For the third term, 

4(cr — l)ua.o + 2\/ZjC(Xq — (4a; — c 2 )ct\ 
>(4o; + 2^/ujc)ol q - (4a; - c 2 )a x 
c 



> 0. 



4w(l H -1=) I /i _ - _1 (cosh(crv / 4o; - c 2 x) - l)dx > 0. 

2v w Jo 
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Thus det[d"(u>, c)] > 0, implying the eigenvalues of d"(u,c) have the 
same sign. Turning to the trace, c 2 — 4(cr — l)u> < c 2 — 4u < for 
a > 2. By (42J), tr[d" (u;, c)] < 0. Hence, the two eigenvalues of d"(co, c) 
are negative. □ 

Closely related to det[d"] is the function 



F{z;a)=(a-1) 



(4.3) 



(cosh y — z) ° dy 



(coshy — z) ° x (z coshy — l)dy 



which helps count the number of positive and negative eigenvalues for 
a G (0,2). Indeed, 

Lemma 4.2. For a G (0,2) and admissible (w,c), det[d"(u, c)] has the 
same sign as F( 



2v^' 



o- . 



Proof. We rewrite (4.1) as 



det[d"(£j,c)] 



2\- — l ---2 



=16((T - l) 2 w 2 ag - [ai(c 2 - 4w) + 2^ca ] 



--16u 2 { (a - iya 2 Q 



2y^ 



cosh(o"\/ 4a; — c 2 x) — \ ) dx 



1 2 



Letting y = ay/Aui — c 2 x, 
(4.4) 

det[d"(w, c)] 



lQu 2 



l a- 2 (l + a)i(Auj-c 2 ) 



-1. ,-i-2 



cr 2 (4o; - c 2 



2v^ 



(7 



□ 



When cr = 1 and 2 G (—1, 1), we have F(z; a) = — 1 and det[d"(a;, c)] = 
— 1/oj. For a G (0, 1) U (1, 2), we can evaluate the function F(z\ a) nu- 
merically, as shown in Figures [l] and [2j For any fixed cr G (1,2), -F(z; cr) 
is monotonically increasing in z and has exactly one root Zq in the in- 
terval (—1, 1). For fixed cr G (0, 1), F(z] cr) is monotonically decreasing 
in z and strictly negative. It is this numerical computation of F which 



is used to complete the proofs of Theorems L3 and L4 In contrast, for 
a > 2, we can prove that F(z; a) is strictly positive without resorting 
to computation. 

We thus have the following theorem about p(d"(u, c)): 

Theorem 4.3 (Numerical). For admissible (u>,c), 



GENERALIZED DNLS SOLITON STABILITY 




Figure 1. (a) Function F(z\ a), Q4.3 ), for several values 
of a G [1,2]. (b) is a magnified plot near the z-axis. 
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= 0.8 
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(b) 



Figure 2. (a) Function F(z; a), ( |4.3 ), for several values 
of a G (0, 1]. (b) is a magnified plot near the z-axis. 
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(i) when a G (1, 2) and c = 2zq^[uj, det[d"(u, c)] = 0, 

(ii) when a G (1,2) andc < 2z§^[uj, det[d"(u,c)} < 0; p(d"(cu, c)) = 1, 
(Hi) when a G (1, 2) and c > 2z y/u, det[d"(u;, c)} > 0; p(d"(u, c)) = 

or p(d"(u, c)) = 2. 
fivj w/ien a G (0, l),det[d" (w, c)} < 0; p(d"(u,c)) = 1. 
fwj w/ien a = 1, det[d"(a;, c)] = < 0; p(d"(u,c)) = 1. 



5. Orbital Stability and Instability 
In this section, we complete the stability/instability proofs. 



Proof of Theorem 1.2 , From Theorem 3.1 and Theorem 4.1, n(H) = 1 
for any a > 0, and p(d") = for a > 2. Thus n(H) —p(d") = 1, is odd 



and all solitary waves are orbitally unstable by Theorem |1.5 



□ 



Proof of Theorem By assumption, for each a G (1,2), there exists 
Zo, a unique zero crossing of ( 4.3[ ), and this function is monotonically 
increasing. By Theorem 4.3, det[d"(u,c)} < for admissible (cj, c) 
satisfying c < 2z y/co. It follows that there is one positive and one 
negative eigenvalues for d"(u,c). Therefore, p(d") = 1. Furthermore, 
from the theorem 3.1, n(H) = 1. Hence 



n(if) -p{d") = 0, 

and we have the orbital stability of solitary waves. 

When c > 2z y/u, also by Theorem 4.3, det[d"(u, c)] > 0. So the 
signs of the two eigenvalues of d"(u,c) are the same. If both of the 
eigenvalues were positive, then 2 = p(d") > n{H^) = 1. This contra- 
dicts (1.17). Hence both of the eigenvalues are negative and p(d") = 0. 



It follows that 

n(H)-p(d") = l, 
and we have the orbital instability of solitary waves. 

Following the same argument, we can prove Theorem |1.4 



□ 



Proof of Theorem l.Jf. , When a G (0, 1], from Theorem 4.3, det[d"(a;, c)] < 
for admissible (uj,c). Consequently, d"(u,c) has one positive eigen- 
value and one negative eigenvalue; p(d") = 1. By Theorem 3.1, n(H) = 
1. Hence 

n(H) -p(d") = 0, 
and the solitary waves are orbital stable. □ 
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6. Discussion and Numerical Illustration 

We have explored the stability and instability of solitons for a gener- 
alized derivative nonlinear Schrodinger equation. We have found that 
for a > 2, all solitons are orbital unstable. Using a numerical calcula- 



tion of the function F(z;a) defined in (4.3). we have also shown that 
for < a < 1, all solitons are orbital stable. For 1 < a < 2, our 
computation of F(z; er) indicates there exist both stable and unstable 
solitons, depending on the values of u and c. In particular, for fixed uj 
and a > 1, there are always both stable and unstable solitons for prop- 
erly selected c. For a near 1, the unstable solitons are always rightward 
moving, but, as Figure [T] shows, the root, Zq, becomes negative as o 
approaches 2. Once zq < 0, unstable solitons can be both rightward 
and leftward moving. 

Other dispersive PDEs possessing both stable and unstable solitons, 
such as NLS and KdV with saturating nonlinearities, [4|[5| [24|[28|[3l] , 
achieve this by introducing a nonlinearity that breaks scaling. In con- 
trast, gDNLS always has a scaling symmetry, and throughout the 
regime 1 < a < 2, the scaling is L 2 supercritical. This also im- 
plies the existence of an entire manifold of critical solitons, precisely 
when c = 2zq\Juj. Along this curve, the standard stability results 
of (7l[8j 34 35], break down, and a more detailed analysis is required. 
In p], the stability can be demonstrated in this degenerate case pro- 
vided d(u, c) remains convex. Given that within any neighborhood of 
a critical soliton there exist unstable solitons, we conjecture that it is 
unstable. While there has been recent work on critical one parameter 



solitons for NLS type equations, |4j[5 24 28 , to the best our knowledge, 
there has not been an analogous work on two parameter solitons. 

While the equation retains the scaling symmetry, we observe that, 
in contrast to NLS solitons, not all gDNLS solitons can be obtained 
from scaling. Indeed, for (1.16), all solitons e lA *i?(x; A), solving 



-AR + XR- \R\ 2a R = 

can be obtained from the A = 1 soliton via the transformation 

e iA *i?(x;A) =e iXt \^R(\^;l). 

In contrast, while the gDNLS solitons also inherit the scaling symmetry 
of gDNLS, not all admissible (u, c) can be scaled to a particular soliton. 
Instead, 

(6.1) ^ lC (x,t) = e*"V w>c (a: - ct) = e^Vi.c/^v 7 ^ - ct)). 
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Only solitons for which 

Cl c 2 



can be scaled into one another. 

Our results were based on the assumption that a weak solution ex- 
isted. While we do not have an H 1 theory in general, our results can, 
in part, be made rigorous as follows. For a > 2, one should be able to 
apply the technique of [33] to obtain a local solution in H s , with s > 1. 



Alternatively, for a > 2 and integer valued, 19,20 can be invoked. 
Again, this yields a local solution in H s , s > 1. For s sufficiently large, 
the solution will also conserve the invariants. 

This is sufficient to fully justify the instability of the unstable soli- 
tons, since there is sufficient regularity such that if the solution leaves a 
neighborhood of the soliton in H 1 , it also leaves in H s , s > 1. However, 
this is insufficient to prove stability, because even if the solution stays 
close in the H 1 norm, the norm of the solution could grow in a higher 
Sobolev norm. 

There is also the question of the monotonicity of F, for which we 
relied on numerical computation for o < 2. Looking at Figures [T] and 
|2j it would appear that the F(z; a = 2) is an upper bound on F(z; a) 
for 1 < a < 2. In addition, there appears to be a singularity at z = 1. 
Likewise, the line F = appears to be an upper bound in the range 
a < 1. A more subtle analysis may permit a rigorous justification of 
our work in this regime. 

Lastly, we provide some numerical experiments of solitons on both 
the stable and unstable branches. We studied the stability near the 
turning point c = Izq^/Zj. When o = 1.5, zq = 0.0618303. The initial 
condition is chosen as 



(6.2) ip (x, 0) = ipuj )C {x, 0) + O.OOOle 



-2x 2 



We simulate ([1/7 ) using the fourth order exponential time difference 
scheme of |17|, and treat the nonlinearity pseudospectrally. Though 
the nonlinearity is not polynomial in its arguments, ip, i/j and ip x , 

H 3 ^x = iH> M fa, 



we found that dealiasing as though it were a quintic problem proved 
robust. 

Our results follows: 



(1) When u = 1 and c = < 2z = 0.1236606, Figure [3j shows that 
the solitary wave retains its shape for a long time (t = 100). 
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^ 1 




100 



Figure 3. Evolution of a perturbed orbitally stable soli- 



ton, u = 1 and c = 0, with initial condition (6.2). 



(2) When u = 1 and c = 0.2 > 1z = 0.1236606, Figure g] shows 
that the amplitude of the solitary wave increases rapidly near 
t = 10 and it is not orbitally stable. 

Our simulation of the unstable soliton suggests that, rather than dis- 
perse or converge to a stable soliton, gDNLS may result in a finite time 
singularity. We will explore the potential for singularity formation in 



the forthcoming work 23 . 



Appendix A. Auxiliary Calculations 

In this section, we present certain integral relations that are helpful 
in studying the determinant and trace of d"(u, c). In the following, we 
denote k = y/Aco — c 2 > and rewrite the solitary solution (1.9) as 

(p(x) 2a = f(cu, c)h(u), c; x) _1 ,with 



/(w,c) 



(a + 1)k 2 
2v^ 



h(x; a; to, c) = cosh(afi:a;) — 



2Ju' 
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Figure 4. Evolution of a perturbed orbitally unstable 



soliton, to = 1 and c = .2, with initial condition (6.2). 



We also rewrite the functionals Q, P defined in (2.6)-(2.7) and their 
derivatives in terms of h and /: 



(A.l) 



(A.2) 



Q 
P 



\p\ 2 dx = f <* i h °dx, 



POD 1 POD 

- / ( p 2 dx + — — / V 2 ° +2 dx 
■ — f° I h~*dx-\ f~^~ I h~^dx. 



(A.3) 
(A.4) 



d c Q =-f ' a ° fc [ h °dx--f° I h ^rh c dx, 



a 



d u Q = -f 1 "" fw I h °dx--f° I h "*'h w dx, 



a 



a 



1 r<x> 
1 „ 1 / . er+1 



a 



2(» 
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1 



C „i- 



d c P = --f- I h~dx-—f~f t 



h °dx 



(A.5) 



c „ i 



+ —f a I h ° h c dx + —f"f c / h " dx 



2a' 

1 l + <7 

f a 

2a 

C „l-o 



2a' 



h + " h r dx 



(A.6) 
where 

fc 



2a' 



c . i 



/ CT fui / h "dx+ — f / h ct h u dx 



2a' 



1+2 



+ -f;L h~~dx - 



c(l + a) _ (1 + a)(4u + c 2 ) 

Jul 



h + " h^dx, 



i CrC • i / \ 1—1/2 7 

h c = x smhlaKX) oo 1 , ti w 

K 2 



2a c _s 

— isinn cr/tx H — a; 2 . 
« 4 



The expressions in (A.3)-(A.6) involve various integrals. The next lem- 
mas show that all of them can be expressed simply in terms of 



dir. 



h—°- n dx. 



First, we have 
Lemma A.l. 

(A.7) 



a 2 



Auj 2cy/u(2 + a) 
-a H t — ; - s - an 



K a + l)n 2 (a + l)n 2 

Proof. We first rewrite «o, and then integrate by parts: 

/■OO ^ poo ^ 

/ h~°~ 1 hdx = — / h~*~ (sxnh.(aKx))'dx — oti 
Jo ° K Jo 



a 



2x/u 



a A 

(CT+I) 



h ° 2 (sin\i 2 (a nx))dx — cti 



2Jco 



a'- 



""" ( " I+ 2^ )2 - 1 ^- ai 2^- 



Regrouping the terms in this last expression, we obtain (A.7). □ 
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Lemma A. 2. We have the following relations: 



(A.8) 
(A.9) 
(A.10) 
(A.ll) 



*h c dx 



CO 



a 2 



POO 

/ h~°~ 1 h c dx 
Jo 



h ° 2 h,,,dx 



h " x h,Ax 



2^j (o + l)n 2 

1 CO 

c 



Qti, 



r"2 + 



2o 



4w 3 /2 z (o + 1)k 2 
c 2o 



Oil, 



Proof. By integration by parts, and n integer, 



h " n h c dx 



o 



i-n+l 



dx- 



Jo 



l dx, 



Choosing n — 2,1, we get (|A.8|) and (|A.9|). The relations (|A.10) and 

i 

□ 



(A.ll) are obtained from (A.8), (A.9) and h, 



— hr — 



K 



4u 3 / 2 



c 



Using Lemmas A.l and |A.2 , we have: 

Lemma A. 3. Denoting k = 2~v~ 2 o~ 1 (l + o)& k 2 ^~ 1 ^oj~^~^ } we have 

d c Q = 2k [2c{o - 2)u 1/2 a + n 2 ai] 
d w Q = kuj- 1 [(2c 2 - 8(<r - u 1/2 a - k 2 c«i] 
d c P = k [(2c 2 - 8(cr - l)u) co l/2 a - ^ca^ 
d^P = 2k [2c(o - 2)cu 1/2 a + K 2 a x ] . 



This is used in Section 4 to obtain (4.1). 
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